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Modeling Neutral Densities Downstream of a Gridded Ion Thruster 


George C. Soulas 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

The details of a model for determining the neutral density downstream of a gridded ion thruster are presented. An 
investigation of the possible sources of neutrals emanating from and surrounding a NEXT ion thruster determined that 
the most significant contributors to the downstream neutral density include discharge chamber neutrals escaping 
through the perforated grids, neutrals escaping from the neutralizer, and vacuum facility background neutrals. For the 
neutral flux through the grids, near- and far-field equations are presented for rigorously determining the neutral density 
downstream of a cylindrical aperture. These equations are integrated into a spherically-domed convex grid geometry 
with a hexagonal array of apertures for determining neutral densities downstream of the ion thruster grids. The neutrals 
escaping from an off-center neutralizer are also modeled assuming diffuse neutral emission from the neutralizer keeper 
orifice. Finally, the effect of the surrounding vacuum facility neutrals is included and assumed to be constant. The 
model is used to predict the neutral density downstream of a NEXT ion thruster with and without neutralizer flow and a 
vacuum facility background pressure. The impacts of past simplifying assumptions for predicting downstream neutral 
densities are also examined for a NEXT ion thruster. 
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Nomenclature 

y-intercept in the non-dimensional wall flux equation 

reservoir particle emission area 

cylindrical wall particle emission area 

slope in the non-dimensional wall flux equation 

neutral particle mean thermal speed 

normalized neutral density distribution upstream of the thruster grids 
charge of an electron (i.e., 1.602xl(T 19 C), 

Maxwell-Boltzmann speed distribution function 
accelerator grid open area fraction 

neutral density correction factor to account for the screen grid (i.e., 0.93 for the NEXT grids) 
thruster beam current 

discharge chamber total doubly-to-singly-charged ion current ratio 
hexagon side number 

Boltzmann’s constant (i.e., 1.3804xl(T 23 J/K) 

transmission probability of neutrals through an accelerator grid aperture 

length between the surface normal of dA 0 to the downstream neutral density location 

grid aperture center-to-center spacing 

length between the grid neutral emission site to the downstream neutral density location 

length between the neutralizer neutral emission area to the downstream neutral density location 

number of apertures along the side of a hexagonal array of holes 

total flow rate into the discharge chamber 

molecular mass of the neutral gas particle 

neutralizer mass flow rate 

downstream neutral density 

far-field equation for the downstream neutral density from an aperture 

far-field equation for the downstream neutral density from a grid by summing individual aperture 
contributions 

far-field equation for the downstream neutral density from a grid by integration of the grid surface 
downstream neutral density from discharge chamber neutrals escaping through the grids 
equation for the downstream neutral density from an aperture in the near- and far-fields 
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nnear-gridsum equation for the downstream neutral density from a grid by summing individual aperture contributions 
that is valid in the near- and far-fields 
n„eut downstream neutral density from neutralizer neutrals 
n 0 reservoir or discharge chamber neutral density 

n o,neut neutralizer reservoir neutral density 

n tank vacuum facility background neutral density 

n tot total downstream neutral density, including contributions from n grid , n neut , and n tank 

n up downstream neutral density from reservoir neutrals through an aperture 

n wa n downstream neutral density reflected from an aperture cylindrical wall in the near- and far-fields 

p number of successive rows defining a hexagonal array of apertures 

p max maximum number of rows defining a hexagonal array of apertures 

Ptank vacuum facility background pressure 

NASA National Aeronautics and Space Administration 

NEXT NASA’s Evolutionary Xenon Thruster 

NSTAR NASA’s Solar Electric Propulsion Technology Applications Readiness program 
R accelerator grid aperture radius 

Rgrid spherically-domed grid radius of curvature 

r downstream neutral density radial location 

r near radial aperture-to-grid coordinate system transformation equation 

r n0 neutralizer keeper orifice radial location 

R no neutralizer keeper orifice radius 

r n0 ff radial offset of the neutralizer keeper emission area 

r Q reservoir radial location 

r 0 f final reservoir, or discharge chamber, radial location limit of integration 

r Q j initial reservoir, or discharge chamber, radial location limit of integration 

r oII radial limit of integration for the reservoir’s, or discharge chamber’s, contribution to the downstream 

neutral density 

r oln radial limit of integration for the reservoir’s, or discharge chamber’s, contribution to the downstream 

neutral density 

t accelerator grid aperture thickness 

T far-field angular distribution function for particles downstream of a cylindrical aperture 

T 0 reservoir or discharge chamber neutral gas temperature 

T tank vacuum facility background particle temperature 

v neutral particle speed 

z downstream neutral density axial location 

Znear axial aperture-to-grid coordinate system transformation equation 

z noff axial offset of the neutralizer emission area 

z w aperture wall axial location 

z wf final aperture wall axial location for integration 

z w ; initial aperture wall axial location for integration 

z wII axial limit of integration for the wall’s contribution to the downstream neutral density 

a angle between the dA 0 surface normal and the downstream neutral density 

oigrid angle between the grid surface normal and the downstream neutral density 

ot n eut angle between neutralizer emission area surface normal and the downstream neutral density 
P downstream azimuthal angle 

y angle along spherically-domed grid 

y max angle along spherically-domed grid to the edge of the perforated grid 

T 0 reservoir or discharge chamber flux 

T w aperture wall flux 

0 neutralizer keeper orifice azimuthal location 

k aperture wall flux variable 

9 cylindrical aperture azimuthal angle 

(p f final cylindrical aperture azimuthal angle limit of integration 

9 i initial cylindrical aperture azimuthal angle limit of integration 
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(Poi angular limit of integration for the reservoir’s, or discharge chamber’s, contribution to the 

downstream neutral density 

9omax angular limit of integration for the reservoir’s, or discharge chamber’s, contribution to the 
downstream neutral density 

<p w i angular limit of integration for the wall’s contribution to the downstream neutral density 

cp wt angular limit of integration for the wall’s contribution to the downstream neutral density 

cq 0 solid angle from dA 0 

co w solid angle from dA w 

Subscripts and Superscripts 

non-dimensionalized value (dimensions are divided by R, densities by n 0 , and fluxes by r o ) 
o reservoir, or discharge chamber 

w aperture wall 

I. Introduction 

The successful demonstration of the NSTAR (NASA’s Solar Electric Propulsion Technology Applications 
Readiness program) ion thruster on the Deep Space 1 mission and its subsequent use on the Dawn spacecraft have 
demonstrated the utility and reliability of ion thrusters for NASA’s deep space missions (Refs. 1 to 5). As a direct 
result of this success, NASA continued the development of ion propulsion technology by investing in NASA’s 
Evolutionary Xenon Thruster, or NEXT, which represents a next-generation ion propulsion system (Ref. 6). The 
objective of the NEXT project is to advance this next generation ion propulsion technology to NASA Technology 
Readiness Level 5, with significant progress towards Level 6, to support NASA Science Mission Directorate 
missions (Ref. 7). The NEXT program is developing an advanced ion propulsion system, that includes, in part, a 
high performance, 7 kW ion thruster with a spherically-domed beam extraction grid system with a chord diameter of 
36 cm. 

Although the NEXT ion thruster exhibits high propellant utilization efficiencies throughout its throttling range, 
some un-ionized propellant, or neutrals, will inevitably escape the ion thruster’s discharge chamber and neutralizer. 
For example, some neutrals will escape from the thruster discharge chamber through the grid apertures and expand 
into the vacuum of space. Some of these neutrals can undergo charge-exchange with the beam ion plume to create 
slow-moving charge-exchange ions. These charge-exchange ions are subsequently directed by potentials within the 
ion beam and charge-exchange plasmas away from the thrust axis and possibly towards spacecraft surfaces. 

Although the energies of these ions are low, modeling of the charge-exchange plasma is typically conducted to 
understand where these charge-exchange ions are directed and to assess their possible impact on the spacecraft. 
Properly modeling charge-exchange production requires an accurate understanding of the neutral density 
distribution downstream of and surrounding an ion thruster. 

In addition to charge-exchange plasma effects on the spacecraft, the neutrals escaping from the thruster can also 
affect plasma diagnostics’ measurements. For example, measurements of the doubly-charged beam ion current 
relative to the singly-charge current are affected by differences in charge-exchange cross-sections (Ref. 8). 
Specifically, the charge-exchange cross-section for singly-charged ions is much larger than that of doubly-charged 
ions (Ref. 9). This can lead to exaggerated doubly-to-singly-charged ion current ratios when measured far 
downstream of the ion thruster and at high background pressures (Ref. 8). Properly modeling this charge-exchange 
effect again requires an accurate understanding of the neutral density distribution downstream of the ion thruster. 

Many past studies have modeled the downstream neutral density emanating from an ion thruster. Unfortunately, 
most studies made simplifying assumptions that limited the utility of their models. These simplifying assumptions 
ignored the effect of the grid shape and the neutral efflux ‘‘‘beaming” effect of individual apertures of finite 
thickness, where thicker aperture walls direct more escaping neutrals along the aperture centerline (Refs. 10 and 11). 
Very early studies either treated the thruster as a point source for far- field analyses or merely assumed that the 
thruster was a flat surface with neutral particle emission throughout the surface, allowing for a simpler analytical 
expression to be used (Refs. 12 and 13). A large number of studies modeled the downstream neutral density using an 
analytical expression that merely assumed a neutrals point source upstream of the thruster grids (Refs. 14 to 21). 
Another study used an analytical expression that assumed that neutrals were emitted from a flat disk diffusely 
(Ref. 22). And yet other analyses utilized a Monte Carlo Particle-in-Cell technique to model the neutral efflux from 
ion thruster (Refs. 23 to 26). Such a technique has the benefit of including the impacts of charge exchange on the 
velocity distribution of the neutral efflux, for example, but computational difficulties with this technique led to 
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simplifications that included computational boundaries that assumed a flat grid without individual apertures. The 
failure of all aforementioned studies to capture the effect of the grid shape and the “beaming” effect of individual 
apertures reduced the accuracy of the modeled results and limited the utility of the models to the far-field. 

Three studies did include the effect of the grid shape and two studies included neutral efflux from individual 
apertures (Refs. 27 to 29). Unfortunately, two studies did not accurately model the neutral efflux “beaming” effect 
of individual apertures (Refs. 27 and 28). The other study correctly captured this neutral “beaming” effect (Ref. 29). 
However, all three studies did not present the equations necessary to recreate the models. 

In addition to these simplifying assumptions, almost all of the aforementioned studies assumed that neutrals 
were only emitted through the thruster grids (Refs. 12 to 27 and 29). So, most of the studies neglected the 
contribution of the neutral efflux from the thruster’s neutralizer and the possible contribution of downstream 
reflections from other thruster surfaces. 

This paper will present the details of a model for determining the neutral density downstream of a gridded ion 
thruster. The model will include the effects of the grid’s spherical shape and the neutral efflux “beaming” through 
cylindrical apertures in both the near- and far-fields. The first section will examine all sources for neutrals 
emanating from an ion thruster. The following section will present near- and far-field equations for modeling neutral 
densities emanating from a discharge chamber through cylindrical apertures based on equations derived in other 
studies (Refs. 10, 1 1, and 30). These equations will then be integrated into the complex perforated grid geometry 
and this, along with the neutral efflux from the neutralizer, will be used to model neutral densities downstream of an 
ion thruster. The model will then be used to predict the neutral density downstream of a NEXT ion thruster with and 
without neutralizer flow and a vacuum facility background pressure. Finally, the impact of past simplifying 
assumptions for predicting downstream neutral densities will be assessed. 

II. Sources of Neutrals From and Surrounding an Ion Thruster 

Figure 1 shows sketches of a typical ion thruster with the various sources of neutrals that can emanate from and 
surround the thruster. As shown in the sketch, the sources of neutrals downstream of an ion thruster originate from 
the following thruster regions: 

1 . un-ionized propellant escaping from the discharge chamber through: 

a. the apertures of the perforated grids; 

b. the thruster’s plasma screen; 

c. the gap between the front mask and accelerator grid; and 

2. un-ionized propellant escaping from the neutralizer; 

3. vacuum facility background neutrals. 




Figure 1. — Unionize propellant escape paths from the thruster (right sketch) and downstream reflections of 
un-ionized propellant escaping from through the grids (left sketch). 
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Also important are the reflections that un-ionized neutrals make after they exit the thruster, as shown in 
Figure 1. It will be shown that the discharge neutrals escaping through the grids are the most significant source of 
neutrals from the thruster. The downstream reflections of these neutrals include: 

1 . neutrals reflected from the front mask; 

2. neutrals reflected from the neutralizer; and 

3. ions that impact the accelerator grid downstream surface and aperture walls and leave as neutrals. 

As the lists above and Figure 1 demonstrate, there are a numerous regions from an ion thruster where neutral 
particles can be ejected. 

Fortunately, there are two considerations that can reduce the complexity of determining downstream neutral 
densities. The first is related to the neutral densities within and surrounding the thruster. Ion thrusters typically 
operate at very low neutral densities. For a NEXT ion thruster operating at full power, for example, the neutral 
density within the discharge chamber is estimated to be about 2.5 x 1 0 1 8 nf 3 . At this low particle density, the mean 
free path for neutral-to-neutral collisions is about 6 cm, which is 1/6 of the NEXT discharge chamber chord diameter 
and about 50x larger than an accelerator aperture diameter. Downstream of the thruster and within a vacuum facility 
at a 6 . 7 x 1 0 3 Pa (5 X 10 6 torr) background pressure, the mean free path for neutral-to-neutral collisions is about 
93 cm. Such large mean free paths define a rarified gas flow regime referred to as free molecular flow. In this flow 
regime, particle-to-particle collisions are so rare that they can be neglected. The benefit of modeling this regime is 
that the resulting downstream neutral density from the various neutral particle sources (shown in Fig. 1) can be 
modeled separately, and the resulting fluxes and densities can be added together to determine the total downstream 
neutral density. 

Only within the neutralizer, where pressures can reach as high as about 67 torr at full power for the NEXT 
neutralizer (Ref. 31), are the mean free paths for neutral-to-neutral collisions small enough for these collisions to be 
significant (i.e., mean free paths of about 3 x 1 0 6 cm at full power within the cathode insert region). However, these 
neutral-to-neutral mean free paths increase rapidly as the gas is ionized (and hence the neutral population decreases) 
and the neutral gas plume expands into vacuum. Because the neutral density downstream of and surrounding the 
thruster is desired, collisional effects within the neutralizer and immediately downstream of the neutralizer cathode 
orifice can also be ignored and the downstream free molecular flow conditions can be considered. 

The second consideration that can reduce the complexity of determining downstream neutral densities is the 
relative contribution of each neutral source and escape path. Although neutrals can escape from many regions of an 
ion thruster, some of the aforementioned emission sites can be neglected without any significant loss of accuracy. 
Although there are three escape paths for neutrals to leave the discharge chamber, the vast majority of particles will 
exit through the grid apertures. This is due to two reasons 1 ) the open area between the grids and discharge chamber 
is only about 1 1 percent that of the accelerator grid open area for a NEXT thruster; and 2) the transmission 
probability through the grid apertures is much higher. Regarding the latter, the escape paths from the discharge 
chamber to the plasma screen and front mask are typically significantly tortuous, involving multiple wall collisions 
before escaping, and exhibiting no direct line-of-sight. However, the escape path through the accelerator grid 
apertures is more direct and, as a result, will have a much higher transmission probability. The transmission 
probabilities through the front mask and plasma screen escape paths were crudely estimated for a NEXT ion 
thruster. This was done by treating the transmission between components within the thruster independently and 
assuming that the entrance conditions of the neutrals for each component was random (Ref. 32). In addition, other 
simplifying assumptions were made to handle bends and elbows, as well as the shapes through which neutrals were 
flowing. The total amount of propellant leaking through the front mask and plasma screen was estimated to be about 
1 percent of that escaping through the accelerator grid apertures. Of this 1 percent neutral loss, most escaped through 
the plasma screen. Because this represents such a small amount of the total leakage, this study will neglect neutrals 
escaping through the front mask gap and plasma screen. 

For a typical NEXT ion thruster, the accelerator current at full power is about 14 mA while the neutral loss rate 
from the discharge chamber is estimated to be 0.51 equivalent Amperes (Ref. 31). So, the total flux rate of 
neutralized accelerator current is only 2.7 percent that of the discharge chamber neutral loss rate. In addition to this, 
the accelerator current collected by the accelerator grid is a combination of charge-exchange ions collected on the 
downstream grid and within the accelerator aperture walls. While those ions collected on the downstream 
accelerator surface would be expected to contribute directly to the downstream neutral density, those collected 
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within the aperture barrel would not necessarily do so. This is because a portion of those ions would be reflected 
upstream and could, therefore, re-enter the discharge chamber. This further reduces the contribution of reflected 
charge-exchange ions on the downstream neutral density. As a result, the total flux rate of neutralized accelerator 
current that contributes to the downstream neutral density is less than 2.7 percent that of the discharge chamber 
neutral loss rate. Because of this small contribution, ions that impact the accelerator grid downstream surface and 
aperture walls and leave as neutrals will also be neglected in this study. 

To determine if neutralizer and front mask reflections are significant, the flux of neutrals that escape from the 
grid apertures and impinge on the neutralizer and front mask was conservatively estimated. It was determined that 
the neutralizer and front mask intercept less than 0.2 percent and less than 0.9 percent, respectively, of the total 
neutral loss from the discharge chamber. Although the reflected neutrals can have a noticeable effect on the neutral 
density very close to the neutralizer and front mask surfaces, the effect further from these surfaces will be negligible. 
The effect of reflected neutrals from the neutralizer and front mask will, therefore, be neglected in this study. 

Vacuum facility background neutrals cannot typically be ignored. Within a vacuum facility at a 6.7 X 1 0 3 Pa 
(5><1 0 6 torr) background pressure, the neutral density would be about 1.6x10 17 nT 3 , which is 6 percent that of the 
NEXT discharge chamber density operating at full power. More importantly, however, the background neutral 
density would be expected to be constant far downstream of an ion thruster, where the density of neutrals emanating 
from a thruster would be expected to be considerably lower than discharge chamber neutral densities. The effect of 
vacuum facility background neutrals will, therefore, be included in this study. 

Neutralizer flow rates are typically large in comparison to the discharge chamber neutral loss rate. For a NEXT 
thruster operating at full power, the neutralizer flow rate is 4.01 seem. Although some of this neutralizer flow will 
exit the neutralizer as ionized gas, the majority is expected to exit as un-ionized propellant. Mikellides modeled a 
NEXT neutralizer during thruster operation at full power (Ref. 33). Fie found that the neutralizer neutral loss rate 
was 2.8 seem during full power operation (Ref. 34). This represents 40 percent of the discharge chamber neutral loss 
rate. Because of the significance of this neutral loss contribution, the neutral loss from the neutralizer will be 
modeled in this study. 

One other source of neutrals, not shown in Figure 1 and neglected in this study, is the fast-moving neutrals 
created by charge-exchange between downstream neutrals and beam ions. Unfortunately, this requires knowledge of 
the spatial distribution of beam ion current density, which is beyond the scope of this study. Because of operation in 
the free molecular flow regime, however, the distribution of these neutrals can be determined by later modeling and 
the result added to that of this study, assuming the change in total neutral density is negligibly small. 

Table 1 summarizes the various neutral escape and reflection sites, their relative contributions, and whether they 
are modeled in this study. As the table shows, the significant contributors to the downstream neutral density include: 
1) the neutral flow from the discharge chamber through the grid apertures; 2) neutrals from the neutralizer; and 
3) background vacuum facility neutrals. The next section will develop the equations for modeling the downstream 
neutral density from these neutral sources. 


TABLE 1.— SUMMARY OF NEUTRAL PARTICLE SOURCES, THEIR RELATIVE CONTRIBUTIONS, 
AND WHETHER THEY ARE MODELED IN THIS STUDY 


Neutral particle source 

Relative contribution 3 

Modeled in this study? 

Discharge chamber — Through grid 
apertures 

99 percent of discharge chamber neutral loss rate 

Yes 

Discharge chamber — Through gap 
between front mask and grid 

~0.2 percent of discharge chamber neutral loss rate 

No 

Discharge chamber — Through plasma 
screen 

-0.8 percent of discharge chamber neutral loss rate 

No 

Neutralizer 

40 percent of discharge chamber neutral loss rate 

Yes 

Facility background neutrals 

6 percent of discharge chamber neutral density 

Yes 

Reflected accelerator current 

2.7 percent of discharge chamber neutral loss rate 

No 

Reflected from neutralizer 

<0.2 percent of discharge chamber neutral loss 
rate reflected 

No 

Reflected from front mask 

<0.9 percent of discharge chamber neutral loss 
rate reflected 

No 

From charge-exchange with beam ions 

Beyond scope of this study 

No 


'Assumes a NEXT ion thruster operating at full power with a background pressure of 5x 10 6 torr. 
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III. Model Development 


The following sections will develop the equations necessary to determine the neutral density downstream of an 
ion thruster. The first section lists the assumptions and fundamental equations used to determine the downstream 
neutral density. Far- and near- field equations for determining the neutral density downstream of a cylindrical 
aperture are presented thereafter. The subsequent sections develop the base equations for determining downstream 
neutral density from a spherically-dome convex grid with a hexagonal array of apertures, an off-center neutralizer, 
and the surrounding vacuum facility neutrals. 

A. Assumptions and Fundamental Equations 

The equations used in the following sections are based on the fundamental equations presented within this 
section. The assumptions used to develop these equations include the following: 

1 . the neutral density is sufficiently low that free molecular flow conditions exist throughout the modeled 
domain; 

2. neutrals from a reservoir (e.g., the discharge chamber) have a Maxwell-Boltzmann speed distribution; 

3. equilibrium conditions exist, so that the rate of neutrals leaving the reservoir is equal to the rate of neutrals 
fed into the reservoir; 

4. all wall reflections are diffuse; 

5. neutral adsorption onto the walls and surface diffusion are negligible; and 

6. wall temperatures are equal to the neutral particle temperature. 

The geometric relation between a differential element from a neutral particle emission site (i.e., a wall or the 
upstream reservoir, such as the discharge chamber) and the downstream neutral density location is shown in 
Figure 2. The variable a is the angle between the neutral emission site surface normal from the differential emitting 
area, dA 0 , to the location of the downstream neutral density, n. The variable L is the distance between the 
differential emitting area and downstream density. 

The downstream neutral density from a reservoir of neutrals can be solved by starting with the Maxwell- 
Boltzmann phase space speed distribution function (Ref. 35). Because equilibrium conditions are assumed, the speed 
distribution function, f(v), is independent of position (i.e., neutral velocities are isotropic). So, the neutral density 
downstream of an orifice to a reservoir is: 


n = J*J*J* n o ' f (v)- d 3 v = J J f(v). n 0 • dra 0 • dv = J f (v)- dv • Jn 0 • dra 0 


0 to. 


Downstream Neutral 
Density, n 


Surface 

Normal 



dA 0 , n 0 

Neutral Emission Site 


Figure 2. — Geometric relation between a neutral 
emission site and the downstream neutral density. 
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Here, n 0 is the neutral density within the reservoir and the differential solid angle, dco 0 , is given by: 


dra 0 


dA 0 -cos(a) 
L 2 


(2) 


The first integral in Equation (1) can be shown to be 1/(4 -ti), so substituting this and Equation (2) into Equation (1) 
yields: 


1 

n = 

4- 71 


f no .£2!fc). dA 

L 


( 3 ) 


This equation is in a form that is useful for solving for the downstream neutral density when the emission site is the 
reservoir. When the emission site is a wall, it is more useful to re-express Equation (3) using a particle flux. This is 
done by noting that the flux in any one direction from a reservoir or diffusely-emitting surface is given by: 



( 4 ) 


where c is the mean thermal speed of the neutrals and the “w” subscript denotes a wall while the “o” subscript 
denotes a reservoir. Solving Equation (4) for n D and substituting this into Equation (1) yields: 


n = 



•f(v)-d 3 v = |f(v)-dv- | 

0 to w 


4-T„ 


•da>„ 


( 5 ) 


Note that T w appears within the second integral because the wall flux will be shown to be a function of the location 
on the wall. Substituting Equation (2) into Equation (5) yields: 

1 f r o \ cos(a) 

n = r w( A w) \ — dA w (< 

c • 71 J L 

A w 

which is in a form that is useful for solving the downstream neutral density when a reflecting wall is the emission 
site. 


B. Equations for the Neutral Density Downstream of a Single Cylindrical Aperture 

The following sections will develop the equations for modeling the neutral density surrounding a grid cylindrical 
aperture. Two sets of equations will be presented. The first is an expression that was initially derived by Clausing 
and is only valid in the far-field (i.e., far enough from an aperture to consider it a point source) (Refs. 10 and 11). 
The second is a series of equations that can be used to rigorously determine neutral density within and downstream 
of a cylindrical aperture both in the near- and far-fields. 

The NEXT ion thruster utilizes a two grid system for accelerating beam ions. The upstream, or screen, grid is a 
thin grid with nearly cylindrical apertures and a high open area fraction. The thin screen grid thickness enhances 
thruster perveance while the small surface area reduces the likelihood that a discharge chamber ion will strike it and 
recombine with an electron, thus improving discharge chamber electrical efficiency. On the other hand, the 
downstream, or accelerator, grid is a thick grid with nearly cylindrical apertures and a much smaller open area 
fraction. The thicker accelerator grid enhances service life capability by providing more material for erosion and the 
smaller open area fraction and increased thickness reduce the neutral transmission probability, thus reducing the 
discharge chamber neutral loss rate. 

The two cylindrical apertures that a neutral particle has to pass through to escape the discharge chamber present 
a difficult problem to solve. Fortunately, though, Kuharski et ah, modeled a two cylindrical aperture set that had a 
geometry that was very similar to that of the NEXT thruster (Ref. 29). They found that the angular distribution of 
particles exiting the accelerator grid was well modeled by assuming a single accelerator grid cylindrical aperture and 
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assuming that the upstream neutral density was merely decreased by 0.93x due to the screen grid. This simplification 
was possible because of the large screen aperture diameter and low thickness relative to that of the accelerator grid. 
This same simplifying assumption will be exploited in this study. 

1. Far-field Equation 

Clausing was the first to develop a far-field expression for the angular distribution of particles downstream of a 
cylindrical aperture (Ref. 10). His equation assumes that the aperture is a point source and that the flux of particles 
from the walls of the aperture is linear as a function of axial thickness, which is a very accurate assumption over a 
range of aperture thickness-to-radius ratios (i.e., t/R) that span from 0 to about 8 (Ref. 30). The angular distribution 
of particles downstream of a cylindrical aperture, T(a), in the far-field is given by (Refs. 10 and 11): 


1- --(1-k)- 


sin p' tan ( a ) | | ^ tan < a ) ^ | t-tan(a) 


1 - 


3-71 


1 f t-tan(a) 

t ■ tan(a) 

2 


if a < arctan 


( 7 ) 


T(a) = 


and 


1-2-k 


3 ■ ji t ■ tan(q) 
2 


if a > arctan| — 
t 


where 


k = - 


4 


t" + 4 - t 


2 + 




(8) 


Here, t is the accelerator aperture thickness and the bar symbol denotes non-dimensionalizing by dividing the 
thickness with the accelerator aperture radius, R. 

Figure 3 shows the geometry for a single aperture that is treated as a point source. It can easily be shown that: 

L = v/r 2 +z 2 (9) 


and 


COS 1 


(a) = 


ir 


(10) 


Equations (9) and (10), along with T(a), can be substituted into Equation (3) to determine the downstream neutral 
density from a cylindrical aperture in the far-field. Note that all variables are independent of dA 0 because a point 
source is assumed, so that integral reduces to jt-R 2 . The resulting equation for the far-field neutral density, n far , is: 


n far( r , Z ) = F S( 


n n R , 

o r r 



( 'Nl 

arccos 

z 

U 2 + z 2 J. 


(r 2 +z 2 ]2 


(11) 


NASA/TM- 


-2010-216930 


9 



Downstream Region 



Upstream Reservoir, 
n 0 


Figure 3. — Geometry for determining downstream 
neutral density in the far-field. 


Here, F scm represents the correction factor to account for the upstream screen grid, which is 0.93 for the NEXT grids 
as discussed earlier. This correction factor can also be set equal to 1 if n 0 is the density upstream of accelerator grid. 
The far-field equation for a cylindrical aperture can be re-expressed in non-dimensional form by: 



The bar symbol above each variable denotes a non-dimensional variable. Neutral density is non-dimensionalized 
with the upstream neutral density, n u , and r and z are non-dimensionalized with the accelerator aperture radius, R. 

2. Near- and Far-field Equations 

Soulas developed equations for rigorously determining neutral density distribution surrounding a cylindrical 
aperture in the free molecular flow regime in the near- and far-fields (Ref. 30). His equations assume that the flux of 
particles from the walls of the aperture is linear as a function of axial thickness. Although he developed equations 
for the upstream, interior, and downstream regions of a cylindrical aperture, only the equations for the interior and 
downstream regions will be presented here. This section will only briefly review the development of these 
equations, although the full equation development can be found in Reference 30. 

The neutral density downstream of a cylindrical aperture is made up of two components: particles that have a 
direct line-of-sight to a downstream location from the reservoir and from the aperture cylindrical walls. In non- 
dimensional form, as indicated by the bar over the variable, the downstream neutral density, n ncar , is given by: 


F rr 

n„ear(rtz) = ^-'< J J — J ' # ' 

r oi 'fi [z 2 + r 2 + — 2 - r - r 0 - cos((j))p 

+ y J rjzj.d-f.cosy 3 # . d - w , 

z»i ii (z - z w ) 2 + r 2 + 1 - 2 • r • cos((|)) 2 


( 13 ) 
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where all dimensions (r, z, t, r Q , and z w ) are non-dimensionalized with the aperture radius, R, and the wall flux, T w , is 
non-dimensionalized with the reservoir flux, r o . The wall flux was assumed to be linear and its non-dimensionalized 
equation is given by (Ref. 30): 


T w (z) = a+b-z 


(14) 


where the above constants a and b are given by: 


a =1 -k 


(15) 


and 


, 2 • k-1 

b = 

t 


(16) 


The first integral on the right hand side of Equation (13) represents the contribution of the reservoir to the 
downstream neutral density while the second integral represents the cylindrical wall’s contribution. The limits of 
integration are a function of the downstream location (i.e., r and z). The contribution of the reservoir and cylindrical 
wall areas depends on which portion of each area has a direct line-of-sight to the downstream location. This is 
illustrated in Figures 4 and 5. As shown in the figures, all reservoir and cylindrical wall areas have a direct line-of- 
sight to downstream locations within region I. But only a portion of the reservoir and cylindrical wall areas have a 
direct line-of-sight to downstream locations within regions II and III. And no reservoir particles have a direct-line of 
sight to downstream locations in region IV. Regions V and IV of Figures 4 and 5, respectively, lie within the walls 
of the cylinder. The reason for the separation of regions II and III of both areas and the further division of 
regions Ilia and Illb for the reservoir was that it was mathematically necessary to correctly integrate these regions 
(Ref. 30). 

Though complicated, the division of the downstream area into the regions illustrated in Figures 4 and 5 were 
used to define a set of equations with different limits of integration that can rigorously determine the neutral density 
downstream of a cylindrical aperture. The base equation for the downstream density, n ncar , in non-dimensional form 
is given by: 


n near 


(r, z) = F scm 



(17) 


where the equations for the upstream (n up ) and wall (n wal i) contributions are listed in Tables 2 and 3, respectively, 
with the appropriate regional border conditions to determine which region, and therefore equation, to use. As before, 
F S cm represents the correction factor to account for the upstream screen grid, which is 0.93 for the NEXT thruster. 
Because of the linear wall flux assumption, the double integrals of Equation (13) could be reduced to single integrals 
in Tables 2 and 3. Note that these equations are applicable to the aperture interior and downstream regions, and that 
the origin for these equations is the upstream entrance to the cylindrical aperture, as shown in Figures 4 and 5, 
unlike the far-field equation. The equations for the variable limits of integration used in Tables 2 and 3 are listed in 
Table 4. 
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Upstream 
Reservoir, n c 




Region I 



Shaded Reservoir 
Portions Have 
Direct Line-of-sight 
to Downstream 
Region 


Figure 4. — Definition of downstream regions for the development of the variable limits of 
integration for the neutral particle flux from an upstream reservoir. 


Upstream 
Reservoir, n 0 



Figure 5. — Definition of downstream regions for the development of the variable limits of 
integration for the neutral particle flux from the cylindrical aperture walls. 
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TABLE 4.— VARIABLE LIMITS OF INTEGRATION USED IN TABLES 2 AND 3 
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C. Convex Grid Geometric Equations and Downstream Neutral Density Equations From a Grid 

The following section will develop the geometric equations for determining downstream neutral density from a 
spherically-domed convex grid with a hexagonal array of apertures. The first part develops the geometric relation 
between the downstream neutral density and the grid surface. The second part develops aperture locations on the 
grid surface for a hexagonal array of holes. A cylindrical coordinate system will be used for the downstream neutral 
density location. 

Figure 6 shows the geometric relation between the downstream neutral density located at (r, z) in a cylindrical 
coordinate system and a neutral emission site located on a spherically-domed convex grid. The angle between the 
normal of the neutral emission surface and the desired downstream neutral density location is given by a gr id and is 
separated by distance L grid . The axial distance, z, is defined from the downstream center of the domed grids. The 
angle P is defined as the angle from the downstream neutral density to the emission site. The variable y is the angle 
from the spherical dome center and y max is the angle to the edge of the perforated grid. The radius of curvature for 
the grid’s spherical dome is given by R gr i d . 

To solve for the downstream neutral density, it is necessary to define a gr i d and L grid as functions ofy, P, R grid , r, 
and z. Using the law of cosines, it can be shown that a gr i d in non-dimensional form is given by: 


a grid ( r , z,y,p) = arccosj 


Z + Rgrid — 

2 

+ Lgrid - (z + Rgrid f ■ tan (y) 2 - f 2 + 2 • (z + Rgrid)- tan(y)- r • cos(P) 

cos(y) gnd 

2- 

Z + R gri d — 

' Lgrid 

cos(y) gnd 


(18) 


where all dimensions are non-dimensionalized with the accelerator aperture radius. Using the law of cosines again, 
the length L grld can be shown to be given in non-dimensional form by: 


"grid 


(r,z,y,p)= yi 


f 2 +Rgnd' 


sin(y) - 2 • r • R grid • sin(y) • cos(p) + {z + R 


'-grid 


• [l-cos(y)]} 2 


(19) 



Figure 6. — The geometric relation between the downstream 
neutral density located at (r, z) and the neutral emission 
site located on spherically-domed convex grid. 
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For both of the aforementioned equations, all dimensions were non-dimensionalized with the accelerator aperture 
radius, R. 

A typical hexagonal pattern of apertures is shown in Figure 7. It will be assumed that a center aperture exists 
within a spherical aperture pattern whose boundary is defined with a maximum arc angle y max (see Fig. 6). All other 
aperture locations can be determined as a function of y and P in spherical coordinates relative to this center aperture 
by adding apertures along the sides of successively larger hexagons centered about this center aperture. Such an 
approach simplifies calculations because of symmetry because the hexagonal pattern can be divided into 6 sides. 

For this study, y and P were defined as functions of three integer values and the aperture center-to-center 
spacing, l cc , explained below and illustrated in Figure 7. Note that for a spherically-domed grid, l cc is an arc length. 
Along a given hexagon’s side, k, and p-l cc distance along a major radius (e.g., along P = 0), where p is the number of 
apertures, aperture centerlines are located m-l cc away from the major radius. It can be shown that y and P as 
functions of k, m, and p are given by: 


P ' ^cc 
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^grid 


• COS 


f 71 

l 6 


i m,p 


COS< 


arctan 


p _ 
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and 
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m,p 
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arctan 
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P 





(21) 



Figure 7. — Typical hexagonal aperture pattern fora spherically- 
domed convex grid. The center aperture is depicted in gray. 
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To determine the y and f:5 locations of every aperture along a hexagon for a given p from the center aperture, y and P 
are calculated from m = 0 to p - 1, and these calculations for P are repeated from k = 1 to 6 to account for all 6 sides 
of the hexagon (because of symmetry, y is independent of k). This process of determining aperture locations in y and 
P coordinates as functions of k, m, and p by building successively larger hexagonal patterns of apertures is continued 
until p-l cc /R gnd < y max . Beyond this, aperture radial locations must be checked against the aperture pattern boundary 
to ensure that apertures at y mp > y max are excluded because they lie beyond this boundary. This process is continued 
until p reaches a maximum value, p max , which is defined as: 


Pmax 


Y max 


J 

R 


cc 

grid 



(22) 


where the brackets about the right hand side of the equation denote the ceiling function. 

With the locations of each individual aperture, a gr jd, and L gri d defined in the above equations, the downstream 
neutral density from a convex spherical dome can now be determined. The far-field neutral density can be solved as 
was done for Equation (12). It can be shown that with the far-field angular distribution function, the neutral density 
downstream of a convex grid, n far _ gridsum , is given in non-dimensional form by: 


n far-gridsum 


(fiz)=- 


F scm ' T [«grid(rt z, y = 0, p = 0)] cos[a grid (r, z,y = 0, p = 0)] 


Pmax 6 p-i 
p-1 k-1 m-0 


4 Lgrid ( r > z,y = 0,p = 0) 

^scm • T [a g ri d (r,z,y m- p,P k m p)] cos[a grid (r,z,y m p ,P k m p )] 

A T 2 “ \~ 11 Ym.p S Yu 

4 L grid (r, z,y m p ,p k mp ) 

and 

^ if Y m.p ^ Y max 


(23) 


The first term on the right hand side of the equation above is the center aperture’s contribution to the downstream 
neutral density while the second term includes all other apertures. The second term on the right hand side also 
includes conditional statements to ensure that apertures located beyond the aperture pattern boundary are not 
included in the calculation. The variables a gr i d , L grid , y, P, and p max are defined in Equations (18) to (22), respectively. 
This equation can be generalized further by accounting for variations in the upstream neutral density within the 
discharge chamber as a function ofy and P (i.e., across the grids). It is assumed that upstream neutral density 
distribution can be described by: 


n o(Y>P) = n o ' d (y,P) 

where D is a distribution function and n 0 is now an average upstream neutral density. So, Equation (23) can be 
generalized with Equation (24) to yield: 
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/— _x _ F scm • T|cc grid (r, z, y = 0, (3 = 0)[ D(y = 0, (3 = 0) cos|a grid (r, z, y = 0, (3 = 0) 

Har-gridsumO’Zj- . ' Tl (- - a n 

Fgrid (b z , Y — P — 0) 
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Pmax 6 p-1 

+m 

p-1 k-1 m-0 


if Y m,p ^ Y 


m,p — i max 


(25) 


and 

0 if Y m,p ^ Y max 


Equations (23) and (25) are non-dimensionalized with the average discharge chamber upstream neutral density, n 0 . 

Further away from the grid, it will be shown that it is reasonable to neglect individual apertures and treat the 
entire grid surface as a neutral-emitting surface. The generalized solution can be solved by starting with 
Equation (3) and noting that: 


dA 0 = R grid -sin(y)-dy -dp 


(26) 


Equations (24) and (26) can be substituted into Equation (3) to yield the neutral density downstream of a convex 
grid, n far _g ridint , in non-dimensional form: 


n far-gridint vE z ) — 
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^grid ' r oa ’ r scm 
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cos[a g r id(r,z,y,P) 
(r,z,y,p) 


T 2 

L grid 
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(27) 


Here, the equation is non-dimensionalized with the average discharge chamber upstream neutral density (n 0 ), a gri d, 
L grid , and T are given by Equations (18), (19), and (7), respectively, and F oa is the accelerator grid open area fraction 
that is given by: 


_ 2-ti R 2 
V3-1 c 2 c 


(28) 


Equation (28) is only valid as long as l cc is very small in comparison to R grid . 

Equation (17), which is valid in the near- and far- fields, can also be applied to the grids. This is done by 
transforming r and z of the aperture coordinate system in Equation (17) to the grid coordinate system (i.e., by 
defining r and z in the aperture coordinate system as functions of y, P, R grid , r, and z in the grid coordinate system of 
Figure 6). The contribution of all apertures can then be summed in a manner similar to Equation (25). It can be 
shown that the result in non-dimensional form is: 


n near-gridsum(rt z ) ^ near [bear (r, z, y = 0, p = 0), z near (r, z, y = 0, p = 0)] ■ D(y = 0, p = 0) 
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(29) 
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where n ncar , y, P, and p max , are given by Equations (17), (20), (21), and (22), respectively, D is given by 
Equation (24), and 


biear ^gridO"’ sinjcXgj^^r, YsP)] (^0) 

and 

^near — ^grid(^i S [cX gritj (r, Z,y,p)]+ t (21) 

The variables a gr id and L gri d are given by Equations (18) and (19), respectively. Equations (29) to (31) are 
non-dimensionalized with the average discharge chamber neutral density, n 0 . 

D. Downstream Neutral Density Equation for a Neutralizer 

For the neutralizer, the gas upstream of the cathode orifice is collisional, so that the aforementioned equations for 
the “beaming” effects of cylindrical apertures in the free molecular flow regime do not apply. Although one option 
is to treat the gas exiting the cathode as a simple choked flow through an orifice, the discharge between the 
neutralizer cathode and its keeper and the ion beam will undoubtedly affect the angular distribution of neutral gas in 
the downstream region. Unfortunately, accurately modeling such effects is beyond the scope of this study. However, 
because the neutralizer discharge would likely randomize the spatial velocity distribution of the neutral gas particles, 
this study will assume that the neutrals exit the neutralizer keeper orifice diffusely. 

Figure 8 shows the geometric relation between the downstream neutral density located at (r, z) in a cylindrical 
coordinate system and the neutralizer keeper orifice. The angle between the normal of the neutralizer keeper orifice 
emission area and the downstream neutral density is given by a neu t and is separated by distance L neut . The axial and 
radial distances, z and r, are referenced to the downstream center of the domed grids, as was done for the grids in 
Figure 6. The angle P is defined as being 0° to the neutralizer centerline. The variables z no ff and r no ff are the axial and 
radial offsets, respectively, of the neutralizer keeper orifice emission area center from the downstream center of the 
domed grid. The variables r no and 0 are the radius and azimuthal angle for the neutralizer keeper orifice area, 
respectively. 

Downstream Location 



Figure 8. — The geometric relation between the downstream neutral density 
located at (r, z) and the neutralizer keeper orifice emission area located on 
the downstream end of the neutralizer. 
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To solve for the downstream neutral density, it is necessary to define a neut and L neut as functions of r, z, P, r no , 
and 0. It can be shown that a neu t is given in non-dimensional form by 


a n eut(r,z,P,r n o,0) = 


arccos- 


z z noff 


L ne ut(r, z ,P,r n o,e) 

The length L neut can be shown to be given in non-dimensional form by: 

Lneut 0"’ Z ’ P> r no ’ ®) — 


(32) 
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r noff r ^ ' r noff ' 


• cos(P) + r n 2 0 - 2 • V r noff + r2 -2-r noff -r-cos(P) • r no • cos(0) + (z - z noff f 


(33) 


For both of the aforementioned equations, all dimensions were non-dimensionalized with the accelerator aperture 
radius, R. 

With a neu t and L neut defined in the above equations, the downstream neutral density from the neutralizer can now 
be determined. The neutral density within the keeper orifice can be solved as was done for Equation (11) and noting 
that: 


n* ncut 


■n-R- 


(34) 


where m neut is the neutralizer mass flow rate, n 0 neut is a reservoir neutral density, R no is the neutralizer keeper orifice 
radius, and m n is the molecular mass of the neutral gas. It can be shown that the neutral density downstream of a 
neutralizer, n neut , is given in non-dimensional form by: 


A neut 


(r,z,p) = 


m n • c • 7T • R 2 , 


2 • 71 R no 

J J 

0 0 


cos[cc neut (r, z, p, r no , 0)] 
Lneut ( r > z > P> r no ■ 


r no • dr no • d0 if z > z 


noff 


and 


(35) 


0 if z < z noff 


The variables a neu t and L neut are given by Equations (32) and (33), respectively. The neutralizer neutral density was 
non-dimensionalized with average discharge chamber upstream neutral density, n 0 . The conditionals ensure that 
there is no contribution from the neutralizer for axial locations less than z no ff because neutralizer neutrals have no 
direct line-of-sight to those axial locations. Note that the neutralizer downstream density is a function of P, unlike 
the equations for the downstream neutral density from the grids, which are nearly axisymmetric so that P can be 
ignored. 

E. Background Vacuum Facility Neutrals 

The background neutrals within the vacuum facility are those particles emitted from the thruster (i.e., neutrals, 
beam ions, charge-exchange ions, and neutralizer ions) that have been reflected from a facility surface. It is assumed 
in this study that free molecular flow conditions exist within the facility and that facility neutral particle velocities 
are isotropic. And although the vacuum facility particles will typically have a lower temperature than those particles 
emanating from the thruster (see, for example, Ref. 36 for NEXT ion thruster temperatures), it will be assumed that 
vacuum facility neutrals reflected from thruster surfaces will still have vacuum facility wall temperatures and are 
reflected diffusely. Finally, it is assumed that the neutral density throughout the vacuum facility is uniform, which is 
a reasonable assumption for thruster operation in large vacuum facilities. 
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Because of the aforementioned assumptions, the inclusion of vacuum facility neutrals is straightforward. The 
background facility neutral density need only be added to the neutral densities from the neutrals escaping from the 
grids and neutralizer. The background facility neutral density, n uank , as a function of pressure and wall temperature is 
given in non-dimensional form by: 


n tank 


n tank 


n 


0 


Ptank 

k B ' Ttank ' n o 


(36) 


Here, p tan k and T mnk are the vacuum facility background pressure and wall temperature and k B is the Boltzmann 
constant. The facility neutral density was non-dimensionalized with average discharge chamber upstream neutral 
density, n„. 


F. Total Neutral Density Downstream of an Ion Thruster 

As discussed earlier, the downstream neutral density contributions of each neutral source can be added together 
because ion thrusters typically operate within the free molecular regime downstream of the thruster. So, the total 
neutral density downstream of an ion thruster, n tot , is given in non-dimensional form by: 

n, 0 , (B z, (3) — Og^lr, z) + n neuk (r, z, P) + n^ an k (37) 

Here, n grid is the discharge chamber contribution given by either Equations (25), (27), or (29); n ncut is the neutralizer 
contribution given by Equation (35); and n tank is the vacuum facility contribution given by Equation (36). Note that 
only the neutralizer density is a function of P, which is 0° in the direction of the neutralizer. Although the hexagonal 
pattern of the apertures within the grid leads to a dependence on p for Equations (25) and (29), it is neglected in this 
equation because the pattern is nearly axisymmetric for the typically small accelerator grid apertures. 

All of the terms within Equation (37) are non-dimensionalized with the average discharge chamber neutral 
density, n 0 . It can be shown that this neutral density is given by: 


1 + 


m,i 


-Ju- 


r; 

2' J + 

T++ 


1 +' 


V '2 ■ 7T • kg • T 0 • m n • RgUd • [l cos(y max )] • F oa • K c • F sc 


(38) 


Here, m d is the total discharge flow rate, Jb is the beam current, e is the charge of an electron (i.e., 1.602x10 19 C), 
J ++ /J + is the discharge chamber total doubly-to-singly-charged ion current ratio, and T 0 is the discharge chamber 
neutral gas temperature. The variable K c is the transmission probability of neutrals through an accelerator grid 
aperture and is given by (Ref. 11): 


n 

I 

K c =2- |T(a)-cos(a)-sin(a)-da (39) 

0 


IV. Results 

The following sections will present the results of the neutral density model developed in this study. Neutral 
densities from a single aperture and a NEXT grid will be presented, and near- and far-field models will be 
compared. The neutral density downstream of a NEXT ion thruster will be presented with and without neutralizer 
flow and a vacuum facility background pressure. Finally, the impact of past simplifying assumptions for predicting 
downstream neutral densities will be assessed. 
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A. Downstream Neutral Density From a Single Aperture 

The neutral density surrounding a single aperture is shown in Figure 9 for a NEXT accelerator grid cylindrical 
aperture. The figure was modeled with Equation (17) and the equations listed in Tables 2 to 4 for the aperture 
interior and downstream regions. Figure 9 also includes additional equations derived in Reference 30 for the 
upstream region. The figure shows a depleted density zone immediately upstream of the aperture, and that the 
neutral density drops rapidly from the upstream region to the downstream region, through the aperture. Figure 10 
shows just the neutral density downstream of the aperture. The neutral density drops to less than 10 percent of the 
upstream value within 1 aperture radius of the downstream surface. 




Figure 10. — Downstream neutral density from a single 
NEXT accelerator gird. 
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The far- field non-dimensional neutral density equation (i.e., Eq. (12)) is compared to that calculated with 
Equation (17) in Figure 11. Here, the neutral density was determined at fixed radial distances from the downstream 
center of the aperture from 2 to 60 radii while the angle was varied (i.e., a spherical coordinate system was used) to 
determine angular distributions. The figure illustrates that at increasing distances from the aperture, the far-field 
angular distribution becomes equal to that of the Equation (17) and that the worst-case error in the far-field equation 
occurs along the axial centerline of the aperture. Also plotted in Figure 1 1 is the angular distribution of a simple 
diffuse assumption (i.e., T(a) = 1) corrected for the reduced transmission probability due to the aperture walls at 
60 radii from the aperture. The diffuse assumption produces a distribution that is more heavily weighted to larger 
angles and results that range from -36 to 107 percent that of the more rigorous calculation. 

Because Equation (17) represents a rigorous solution for downstream neutral density, the far-field equation of 
Equation (12) was compared to it in Figure 12 to determine the error as a function of non-dimensional radial 
distance from the downstream aperture center at various angles for a NEXT accelerator grid aperture. The figure 
shows that maximum far-field equation errors reduce to about 5 percent at about 25 aperture radii downstream for 
the 0° case, but within 9 aperture radii for the remaining angles. This figure demonstrates that the far-field equation 
will produce accurate results within 9 aperture radii downstream of a NEXT accelerator grid aperture. 



Figure 11. — Non-dimensional neutral density angular distributions for a 
variety of radial distances in a spherical coordinate system that is 
centered at the downstream aperture center for a NEXT accelerator 
aperture. Solid lines were determined from Equation (17) and dashed 
lines from Equation (12). 



Figure 12. — Far-field error as a function of radial distance at various 
angles in a spherical coordinate system that is centered at the 
downstream aperture center for a NEXT accelerator aperture. 
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B. Downstream Neutral Density From the Ion Thruster Grids Only 

For this section, only the discharge chamber neutral loss through the grids is modeled. The contributions of the 
neutralizer and background vacuum facility neutrals are neglected so that near- and far-field comparisons can be 
made. In addition, the neutral density within the discharge chamber was assumed to be constant, so that D(y, P) = 1 
in Equation (24), and the effect of the upstream screen grid was included (i.e., F scm = 0.93). 

The non-dimensional neutral density downstream of a NEXT ion thruster from neutrals escaping through the 
grid apertures is shown in Figure 13. This was determined using Equation (29), which includes the rigorous 
calculation of Equation (17) for the downstream neutral density from a single aperture. As the figure illustrates, the 
neutral density drops rapidly from the ion thruster grids into the downstream region due to the rapid expansion of the 
neutrals as they exit the accelerator apertures. 

Close to the thruster grids, the contribution of adjacent apertures to the downstream neutral density of any given 
aperture can be examined. The neutral density downstream of the center aperture of a NEXT accelerator grid (i.e., 
r = 0) from Equation (29) and the adjacent webbing (i.e., r = l cc /2) is compared in Figure 14 to that from a single 
aperture from Equation (17). The number densities at r = 0 and r = l cc /2 converge within 3 aperture radii for the 
multiple aperture case and within 6 radii for the single aperture case. More important, however, are the downstream 
magnitudes. The multiple apertures of the grids increase the downstream neutral density of the center aperture by 
lOx compared to that of the single aperture at only about 6.5 aperture radii from the accelerator grid. This 
demonstrates the significant contribution of adjacent apertures to the downstream neutral density of any given 
aperture for the NEXT ion thruster grids. 

Figure 15 also demonstrates another simplifying assumption that can be employed for the NEXT grids. 

Equation (27) neglects the contribution of individual apertures by treating the entire domed grid surface as a neutral- 
emitting surface whose angular distribution is corrected for the “beaming” effects of the cylindrical apertures. 
Because differences between the densities of Equations (27) and (29) decrease to 1 percent within about 5 
accelerator aperture radii from the accelerator grid, such a simplifying assumption also produces accurate neutral 
densities within about 3 mm of the NEXT thruster’s accelerator grid. Because of the ease of using this equation and 
the accuracy of the results, Equation (27) will be used for the remainder of this study. 



Figure 13. — Non-dimensional neutral density downstream of a NEXT 
ion thruster. The neutral density within the discharge chamber is 
assumed constant. 
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Figure 14. — Non-dimensional downstream neutral density from a NEXT 
ion thruster along the grid centerline and at r = l cc /2 comparing the 
effect of the multiple apertures of a grid to that of a single aperture. 



Figure 15. — Non-dimensional downstream neutral density from a NEXT 
ion thruster along the grid centerline and at r = l cc /2 using 
Equations (25), (27), and (29). 

Although Equation (17) represents a rigorous solution for downstream neutral density of a single aperture, its use 
for detennining the downstream densities from a perforated grid in Equation (29) can be cumbersome and 
computationally costly. Elowever, the far- field assumptions used for solving the downstream neutral density of a 
perforated grid in Equations (25) and (27) are significantly less cumbersome to use and can be computed more quickly. 
So, although Equations (25) and (27) are accurate in the far-field, it would be of value to determine how accurate they 
are close to a NEXT ion thruster. Figure 15 compares the results of Equations (25) and (27) to Equation (29) along the 
centerline of the grids and at a radius of l cc /2 (i.e., between two apertures). As the figure illustrates, the far-field results 
converge rapidly onto the more rigorous results of Equation (29). In fact, differences reduce to 1 percent within about 5 
accelerator aperture radii from the accelerator grid. Because the NEXT accelerator aperture diameter is so small in 
comparison to the overall grid size, the far-field assumptions used in Equations (25) and (27) produce accurate neutral 
densities within about 3 mm of the NEXT thruster’s accelerator grid. 
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C. Total Downstream Neutral Density 

The contributions of the neutralizer and background vacuum facility neutrals are included in this section. For all 
of the following calculations, a NEXT ion thruster geometry was assumed. In addition, the neutral density within the 
discharge chamber was assumed to be constant so that D(y, P) = 1 in Equation (24), and the effect of the upstream 
screen grid was included (i.e., F scrn = 0.93). 

The effect of neutrals from the neutralizer on the downstream neutral density for an NEXT ion thruster operating 
at full power is shown in Figure 16, which excludes the effect of background vacuum facility neutrals and therefore 
represents operation in space. For a NEXT thruster operating at full power, the neutral loss rates are estimated to be 
7.05 seem for the discharge chamber and 2.8 seem for the neutralizer. The discharge chamber neutral temperature 
was assumed to be 300 °C while the neutralizer neutral temperature was assumed to be 1100 °C (Ref. 36). The 
figure illustrates the impact of the neutralizer neutrals on the downstream neutral density. The high neutral density 
surrounding the neutralizer was due to neutrals being distributed across such a small keeper orifice area compared to 
that of the grid area. Figure 17 plots the non-dimensional neutral density as a function of axial location with and 
without neutralizer neutrals at radial locations of 18 cm near and opposite the neutralizer (i.e., P = 0° and 180°, 
respectively). The 18 cm radial location was selected because it is the edge of the perforate grid radius. As expected, 
the neutralizer does not contribute to the density upstream of the neutralizer tip because there is no direct line-of- 
sight there. However, beyond that, the density increases by as much as 68 percent close to the neutralizer. Further 
downstream, the neutralizer’s off-center effect is minimized, but the density is still about 22 percent larger because 
of the neutralizer’s contribution. 

The effect of background neutrals from the vacuum facility on the downstream neutral density for an NEXT ion 
thruster operating at full power is shown in Figure 18. Here, the background pressure was set to 6.7x10 Pa 
(5 X 1 0 6 torr) with a background neutral temperature of 20 °C. The effect of neutralizer neutrals was included in this 
figure using the same assumptions as for Figure 16. The figure illustrates the significant impact of the facility 
neutrals on the downstream neutral density. Figure 17 includes the effect of vacuum facility neutrals on the non- 
dimensional neutral density as a function of axial location at radial locations of 18 cm near and opposite the 
neutralizer. Although the contributions of the grid and neutralizer are significant close to them, the vacuum facility 
contribution represents the vast majority of the neutral density further away. 



Figure 16. — Non-dimensional neutral density downstream of a NEXT ion thruster 
operating at full power in space. The neutral density within the discharge 
chamber was assumed constant and the neutral loss rate from the neutralizer 
was 2.8 seem. 
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Figure 17. — Non-dimensional neutral density downstream of a NEXT 
ion thruster operating at full power with and without neutralizer and 
vacuum facility contributions. The two radial locations shown include 
r = 18 cm in the direction of the neutralizer (i.e. , p = 0°) and r = 18 cm 
opposite the neutralizer (i.e., p = 180°). The neutralizer tip is located 
at about z = 6 cm. 



Figure 1 8. — Non-dimensional neutral density downstream of a 
NEXT ion thruster operating at full power with vacuum facility 
neutrals. The neutral density within the discharge chamber 
was assumed constant and the neutral loss rate from the 
neutralizer was 2.8 seem. Vacuum facility background 
pressure was set to 6.7x1CT 3 Pa (5*1Cr 6 torr). 
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D. Assessment of Past Assumptions 

This section will assess some past simplifying assumptions that were made to more easily determine the 
downstream neutral density. These simplifying assumptions include neglecting the effect of neutralizer neutrals, 
neglecting the “beaming” effect of the cylindrical aperture walls, and assuming a flat grid. All of these assumptions 
will be examined with a NEXT ion thruster operating at full power. 

The impact of the neutralizer was already assessed in the prior section (see Fig. 17). It was shown that 
neglecting the neutralizer’s contribution to the downstream neutral density would lead to a neutral density that is 
22 percent lower than actual in the far-field for a NEXT thruster operating at full power. The impact of neglecting 
the “beaming” effect of the aperture walls is illustrated in Figure 19 along the grid centerline. In this figure, the 
effects of neutralizer and vacuum facility background neutrals were ignored to facilitate direct comparisons. And 
although neutrals were assumed to exit an aperture diffusely, the upstream neutral density was set equal to that 
where the “beaming” effect is included (i.e., the baseline case). The figure shows that the neutral density begins 
higher than the baseline case close to the grids because the diffuse angular distribution contributes more neutrals to 
the grid centerline close to the grid. But further away, this same effect leads to lower neutral densities as the neutrals 
expand more rapidly. So, ignoring the “beaming” effect of the aperture walls leads to an overestimation of 
14 percent close to the grid and an under-estimation that approaches 30 percent in the far- field. 

The impact of assuming a flat grid is also illustrated in Figure 19 along the grid centerline. The upstream neutral 
density was set equal to that of the domed grid case, and the “beaming” effect of the aperture walls is included. The 
figure shows that the neutral density is always higher than the baseline case because the flat grid collimates the 
outgoing neutrals. The flat grid assumption led to neutral density overestimations along the grid centerline of as 
large as 29 percent in the near-field and 14 percent in the far- field. 



Axial Location, cm 

Figure 19. — Non-dimensional downstream neutral density and errors 
as a function of axial location along the grid centerline for: a domed 
NEXT grid that includes the “beaming” effect of the aperture walls 
(Clausing, or the baseline case); a domed NEXT grid that excludes 
the “beaming” effect of the aperture walls (diffuse assumption); and 
a flat grid that includes the “beaming” effect of the aperture walls. 
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Conclusions 


The details of a model for determining the neutral density downstream of a gridded ion thruster were presented. 
A rigorous examination of the possible sources of neutrals emanating from and surrounding a NEXT ion thruster 
concluded that the most significant contributors to the downstream neutral density included discharge chamber 
neutrals escaping through the perforated grids, neutrals escaping from the neutralizer, and vacuum facility 
background neutrals. Because ion thrusters were shown to operate in the free molecular flow regime, each neutral 
density source was modeled separately and the results added together. 

For the neutral flux through the grids, near- and far-field equations were presented for rigorously determining the 
neutral density downstream of a cylindrical aperture. These equations included the “beaming” effect of the aperture 
walls. These equations were subsequently integrated into a spherically-domed convex grid geometry with a 
hexagonal array of apertures for determining neutral densities downstream of the ion thruster grids. The neutrals 
escaping from an off-center neutralizer were also modeled assuming diffuse neutral emission from the neutralizer 
keeper orifice. Finally, the effect of the surrounding vacuum facility neutrals was included and assumed to be 
constant. 

The model was then used to predict the neutral density downstream of a NEXT ion thruster with and without 
neutralizer flow and a vacuum facility background pressure. Individual aperture results showed that far-field 
equation errors reduced to about 5 percent within 9 aperture radii, indicating that the far-field equation will produce 
accurate results close to a NEXT accelerator grid aperture. The neutral density from the NEXT thruster grids as a 
function of axial distance was found to decrease rapidly due to the rapid expansion of the neutrals as they exit the 
accelerator apertures. Because the NEXT accelerator aperture diameter is so small in comparison to the overall grid 
size, the far-field equations were shown to produce accurate neutral densities within about 3 mm of the NEXT 
thruster’s accelerator grid. One of these far-field equations neglected the contribution of individual apertures by 
treating the entire domed grid surface as a neutral-emitting surface whose angular distribution is corrected for the 
“beaming” effects of the cylindrical apertures. Because this equation was least cumbersome and computationally 
efficient, it was used for the remainder of the study. 

The impacts of the neutralizer neutral loss rate and the surrounding vacuum facility neutrals were separately 
assessed for a NEXT ion thruster operating at full power. Beyond the neutralizer, the density increased by as much 
as 68 percent at a radial location of 18 cm close to the neutralizer. Further downstream, the neutralizer’s off-center 
effect was minimized, but the density was still about 22 percent larger because of the neutralizer’s contribution. 
Proper modeling of charge-exchange production in the ion beam requires incorporating the effects of neutralizer 
neutral efflux. Vacuum facility effects were assessed at a background pressure set to 6.7><1 0 3 Pa (5x1 0 6 torr). 
Although the contributions of the grid and neutralizer were still significant close to them, the vacuum facility 
contribution represented the vast majority of the neutral density further away. 

The impact of past simplifying assumptions for predicting downstream neutral densities were also assessed for a 
NEXT ion thruster. Neglecting the “beaming” effect of the aperture walls along the grid centerline led to an 
overestimation of neutral density by 14 percent close to the grid and an under-estimation that approached 30 percent 
in the far-field. Assuming a flat grid led to neutral density overestimations along the grid centerline of as large as 
29 percent in the near-field and 14 percent in the far-field. 
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